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A Three-Dimensional Viscoelastic Model of a Strike Slip Fault 


Summary 


An analytic approximation to the Green's function for the 
displacements due to a strike slip point source in an elastic 
layer over a viscoelastic half-space is developed. This approxi- 
mate Green's f inaction is useful because ir can be analytically 
integrated over the fault surface. Comparison with a numeri- 
cal integration of the exact solution integral indicates that 
the approximation is quite good. The approximate Green's function 
is integrated analytically to obtain the displacements due to 
a finite rectangular strike slip fault in an elastic layer over 
a viscoelastic half-space. Ground displacements and angle changes 
from a model survey net are computed to illustrate the visco- 
elastic relaxation which follows a fracture m the elastic region. 
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Recently there has been a growing interest in the use of 
static or quasistatic displacements , strains and tilts for the 
investigation of earthquake-related phenomena. In particular, 
aseismic horizontal surface displacements due to the San Andreas 
fault have been studied extensively (Savage & Burford 1973; 
Thatcher 1974, 19 75a, b) by the use of geodetic triangulation 
techniques. The importance of the San Andreas lies in its well- 
documented history (Meade 1973) of aseismic surface motions. 

With this data we can learn much about the processes involved 
in pre and postseismic displacements and in tectonic plate motion. 

Observations of surface movements due to the San Andreas 
can be interpreted by means of mathematical models obtained 
by idealizing the fault as a vertical fracture in an elastic 
medium. The models describing near field displacements are 
obtained by integration of a Green’s function for a strike slip 
source. For convenience, the source is usually assumed to be 
located in an elastic half-space whose elastic properties may 
vary with depth. Although the Green's function has been com- 
puted for both line and point sources in layered half-spaces 
(Ben-Menahem & Singh 1968; Rybicki 1971; Chinnery & Jovannovich 
1972; Jovannovich et al. 1974a ,b) the Green's function for a 
point source has nor been analytically integrated over a fault 
surface. The .line source Green's function is really only use- 
ful in describing a fault whose dislocation function changes 
slowly along its strike. The point source Green's function 
can be used to describe a finite-dimensional fault with an arbi- 
trary dislocation function. 
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In this paper we construct a model for a rectangular verti- 
cal strike slip fault which can be used to explain the post- 
seisraic surface displacements occurring after a large earth- 
quake. We propose to represent the fault as a fracture in an 
elastic layer over a viscoelastic half-space, a model conceptually 
similar to that proposed by Nur and Mavko (1974) for infinitely 
long, dipping thrust faults. Unlike the Nur-Mavko model how- 
ever, ours describes a vertical, finite-dimensional strike slip 
fault. In both models, the displacements from a fracture in 
an elastic layer overlying an elastic half-space are first com- 
puted, and the correspondence principle (Biot 1954; Lee 1955; 

Fung 1965) is then used to introduce Newtonian viscoelastic 
properties into the half-space. 

We base our model on an approximate representation of the 
Green' s 'function for horizontal displacements due to a strike 
slip point source in an elastic layer over an elastic half-space. 
Although the exact expression for the Green's function has been 
found (Ben-Menahem & Singh 1968) , it is in the form of an inte- 
gral over a Fourier wave number k. In our model, three additional 
integrations of this Green's function are needed: two over 

the fault surface and an inverse Laplace transform from use 
of the correspondence principle. These four integrals can be 
done numerically, but since near field displacements are desired 
over a wide range of times, small mesh sizes are needed in the inte- 
grations. The approximate Green's function described below 
is accurate, easy to use and has the advantage that all the 
succeeding integrations can be done analytically. Thus the 
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inverse problem (Rundle & Jackson 1976a,b) can be done simply 
using analytically computed partial derivatives. 



An Approximate Image Technique 


Xshii and Takagi (1967) solved a two-dimensional problem 
for the displacements from a line force located at the coordinates 
x 3'' x 2 = an i n fi n ite elastic medium whose rigidity changes 

abruptly at y 3 = 0 (figure 1) . Sybicki (1971) showed that the 
equations of elasticity for this case reduce to Poisson's equa- 
tion. He was thus able to use the method of images (figure 2) 
to find the displacements from a line source in and below a 
horizontal layer over a half-space. For an infinite medium 
whose Y]_ ~ y 2 plane divides the space into two parts, one of 
rigidity ]i ^ and another of rigidity ]i 2 , we can write the dis- 
placements in the y^-direction at a point (y 2 , y^) due to a 
line force at (x 2 , x g) unit magnitude: 
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where is the rigidity for y^ < 0 and p 2 is the rigidity for 
y^> 0. The superscript of >K^ indicates the coordinate along 
which the displacement appears , and the subscript indicates 
the direction along which the force acts. The tilde under the 
variable denotes a vector. 

r 

It has been shown (Rybicki 1971} that the displacements 
V 12 at ^ x 2 7 x 3^ ^ ue to an elementary single-couple line dis- 
location U dS at (y 9 , y_) across a vertical plane can be obtained 
from k:jh;£, by the operation 

v L ( S' v - v s ^ K i ( ^' 

pi y, < o 

v = { 1 J (2) 

**2 y 3 >0 

where is the magnitude of the dislocation and dS is the area 
over which it occurs. It is important- 'to notice that the dis-. 
placements £) at in direction k due to a strike slip 

source at ^ are proportional to the horizontal (12) shear stresses 
at % due to a unit force at x in direction k. 

We would like to find the displacements due to a line dis- 
location in a layered half-space. Using the coordinate sys- 
tem of figure. 2 , one can superpose image line forces to satisfy 
the requirements of a traction-free surface at y^ - 0 and con- 
tinuity of the solution and its consequent stresses at y^ = H. 

The displacements from a line force m a layered half-space 
whose layer thickness is H (Rybicki 1971) is: 
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where is the layer rigidity and" u 2 is the half-space rigidity. 

The operator U dsu^ — • can again be used to obtain the displace- 

0 o y 2 

ments at (x 2 , x 0 ) due to a line strike slip source at (y 2 , y 3 ) . 
Chinnery and Jovannovich (1972) have extended (3) to a many- 
layered half-space. 

Equation (3) satisfies the conditions which uniquely define 
the Green's function for the linear differential operator, the 
Navier equation of elastostatic equilibrium, and the boundary 
conditions of the problem (Courant i Hilbert L953, p. 353) : 

(i) For fixed x 2 and x 3 , G^(£ ; £) has a singularity 
in the stress at y 2 = x 2' Y 3 = X 3* Tlms in two dimensions G^Y, x ) 
is proportional to the logarithm of the distance between observa- 
tion and source points. 

(ii) The displacements G^(^, x) and tractions on any 
horizontal plane are everywhere continuous except at Y 2 = ^ 2 , 

Y3 = x 3 * 

1 

(lii) G^(;£, & satisfies the boundary conditions. 

These three requirements imply that for the Green's function 

1 1 

considered here, the reciprocal theorem holds: G-^ ('/, = G-^(^, ^) . 

Since we are interested in calculating displacements due 
to a finite-dimensional strike slip fault in an elastic layer 
over a viscoelastic half-space we must find the Green's function 
for a strike slip point source. Let us consider the situation 
in which a- double-couple source at y^ is in the layer and attempt 
to write down an approximate Green's f inaction. The’ approxi- 
mation must satisfy conditions (i) - (iii) as well as possible. 

For example, one can show that condition (ii) can be satisfied 
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using any function whose y 3 and x 3 arguments appear only in 

the form {y 3 -x 3 ) 2 / (y 3 + x 3 ) 2 , (y 3 ± 2mH ± x 3 ) 2 , with the function 

occurring in the same term-by-term arrangement as the logarithms 

in equation (3) . This suggests that we can apply the form of 

(3) to our three-dimensional problem. 

We already know the exact analytic solution for the dis- 
1 

placements u^ 2 (;g, Jg) from a double-couple force in an infinite 
medium (Steketee 1958) : 
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The double-couple source is located at (0, 0, y^) and the 
observation point is at U Q dS is the elementary displacement 
dislocation corresponding to the double-couple force, a = ^ , 

g 2 2 

and R = /x^ + x_ 2 + (x 3 - y ) . Again the superscript indicates 
the displacement direction, and the subscripts correspond to 
the fact that the u^ 2 are derived by finding the horizontal 
shear stresses at jr due to a point force at along coordinate 
k. Using the fact that equation (3) represents horizontal dis- 
placements due to horizontal forces, we can try the following 



as an approximate Green’s function: 
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where again is the layer rigidity and y 3 iS the half-space 
rigidity. 

This trial Green’s function has a singularity at the pro- 
per point; its contribution to displacement is continuous at 

3W 12 

x = H; and its contribution to shearing stress, y-g^. 1 snd nor 

5 9W 12 • 1 
mal stress, A-g^ — is continuous at x^ = H if A^ = A ^ , which 

we shall henceforth assume. However, it is evident that there * 
still exists a normal stress at = 0, and thus the surface 
is not stress-free. We can remedy this defect by adding to 
each source and its image a term corresponding to the displace- 
ments that result from superposing a normal stress to cancel 
the normal stress mentioned above: 
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Note that each term of the form — j(A + + — jB + ) is the 

r . r 

Green's function for a strike slip point source at depth 2mH ± y^ 
in a homogeneous half-space for observation points at the free 
surface x 3 = 0 (Steketee 1958) . Thus for observation points 
x 3 > 0, g^ 2 ;£) can be written out by substituting for each 

Green's function for x 3 = 0 the corresponding Green's function 
fo:r arbitrary Xg as found by Steketee (1958) . 

The surface x 3 = 0 is now stress-free, but g^ 2 an ^ 
its associated contribution to the stresses are no longer con- 
tinuous at x 3 — H. We are missing a set of terms which "fixes 
up" 9^2 ^ and makes it smoothly continuous across the boun- 
dary between layer and half-space. Since we are limited to 



the observation of surface motion however, it is necessary only 
to insure that this approximate solution is nearly equal to 
the real solution for observation points at = 0. 

We shall postpone further discussion of the accuracy of 
the approximate Green’s function to a later section. Instead 
we complete the development of the model by integrating (6) 
over a rectangular fault surface and introducing viscoelastic 
properties into the half-space. 
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The Finite Fault 


The solution for horizontal displacements due to a finite 
fault can be found by integrating (6) over the area of the fault 
according to the well-known Volterra relation (Steketee 1958) : 

“ "fault surfaced ?> n e aS <S> (7> 

where n is the outward normal to the fault surface and A u. 

e i 

is the displacement jump. Both a lower tilde and index nota- 
tion are used to denote vectors and the summation convention 
is assumed. An upper carat denotes a unit vector. 

Note that the observation point is. now located at and 
the integral is performed over the variable Performing (7) 

over a rectangular strike slip fault of semilength L with top 

at depth d, bottom at depth D and with Au. = U x. = constant 

i _ o i 

is particularly easy. The result can be taken from the litera- 
ture (Chinnery 1361, 1963) using the coordinate system of figure 
3. For y 3 = 0, x 3 £H, the case of interest here. 
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Introduction of Time Dependence 


Having found the displacements due to a finite-dimensional 
fault m an elastic layer of rigidity ji over an elastic half- 
space of rigidity y ^ , we modify the solution to add Newtonian 
viscous properties to the half-space. It is not our purpose 
here to enter the debate upon the rheological properties of 
the earth’s consitutive material. We remark, however, that 
a Newtonian flow law may be justified for the shallow visco- 
elastic regions in the earth (McKenzie 1968) although there 
is evidence for nonlinear rheology as well (Post & Griggs 1973; 
Stocker & Ashby 1973) . In addition, it has been shown that 
under some conditions a Newtonian rheology produces > f low.istruc- 
tures very similar to a rheology in which the strain rate depends 
upon the cube of the deviatoric stress (Parmentier et al. 1976) . 
Thus we conclude that the use of a Newtonian flow rheology for 
our half-space is justified provided it is recognized that the 
viscosity may be an "average" in the sense of Parmentier et 
al. (1976) . 

It may also be that mineral rheology plays little role 
in the stress relaxation process. Since the zone of brittle 
fracture, ends at. a relatively shallow depth on the San Andreas 
(Brace & Byerlee 1970) , an alternative source of the viscous 
properties of the lower viscoelastic zone may be transient flow 
of water in porous rock (e.g. Nur & Schultz 1973). Darcy’s 
law, which governs hydrogeologic flow, implies that the stress 
relaxation process is essentially linear. Thus to a good 
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approximation, the half-space can be assumed to be a linear 
viscoelastic solid. 

To add viscous properties to the half-space, we use the 
correspondence principle (Biot 1954; Lee 1955; Fung 1965} as 
outlined in the introduction. This principle states that if 
the elastic solution to a problem is known and the inertial 
forces are negligible, the quasistatic solution for a linear 
viscoelastic medium is obtained by replacing all time-dependent 
quantities by their Laplace transforms and then taking the inverse 
transform of the resulting expression. In the case of (6) or 
(8) , we replace p 2 by u 2 (s) and 0 o H(t) by U Q /s, where H(t) 
is the Heaviside step function, t is time and s is the conju- 
gate Laplace variable. 

We choose the properties of the half-space to be those 
of a Maxwell viscoelastic solid whose elastic constants are 
the same as those of the layer. This choice of Maxwell pro- 
perties is motivated by the observation that for times short 
compared to a year rhe earth behaves elastically, while for 
longer times permanent non-elastic deformations are observed, 
for example in. areas of crustal rebound (McConnell 1965; Crittenden 
1968. peltxer, i.974) and in tectonically active fault zones 
during recent geologic history. 

The Maxwell solid is the simplest kind of medium to exhibit 
this dual behavior. A Maxwell element is shown in figure 4; 
if a constant displacement is suddenly applied to the free end, 
the spring will immediately stretch. In the course of time, 
the spring will gradually unstretch as the piston in the dashpot 
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moves . The constitutive equation for a Maxwell solid is 
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where a. . is the stress tensor, e . is the srrain tensor, §. . 
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is the Kronecker delta, u q is a constant rigidity and g is the 
viscosity. The constants A. and y are related to^compres— 
sional stresses which are assumed not to relax. 

Upon performing the Laplace transform on (10) , we can set 
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x is a characteristic time constant for the system. Since the 
rigidity of the Maxwell medium has been chosen equal to that 
of the layer, only the effective viscosity g remains unknown. 

After insertion of (11) in (8) and (9) we must find the 
inverse Laplace- transform l" 1 ^ ^ of quantities such as 
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From Abramowitz and Stegun (1970) we find 
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Upon insertion of (13) into (8) and (9) we get expressions 
of the form 
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where the (^) and B jn (^) can be found from (8) and (9) and 

the F. (7/) are the uniform elastic half-space solutions. 
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Upon examination of (14) it is clear that the first response 
of the system is entirely elastic. As time proceeds,- succes- 
sive terms under the summation on m become progressively more 
important. In the limit t-* 00 , the exponential dies away and 
we get 
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Numerical tests indicate that only the first five terms need 
be retained under the sum in (14) for convergence accuracy within 
1% for tjc2x. Since A m+ (^j and (y) have the same sign as 
F^(y) respectively F^ (^) , we arrive at the conclusion that with 
the- viscoelascic model (14) the quasistatic stress relaxation 
tends to increase the displacements found immediately after 
the fracture. 

For a long fault which breaks entirely through the elastic 



layer, the post-fracture viscoelastic displacements can be an 
appreciable fraction of the initial movement {figures 5-9). 

Within ten kilometers of the fault shown, there is a relatively 
small amount of post-fracture movement compared to the rupture- 
induced displacements. Between approximately ten and one hun- 
dred kilometers, the post-fracture displacements diminish slowly 
with distance at any given time. At distances of more than 
thirty kilometers from the fault, most of the total movement 
is accounted for by viscoelastic effects rather than by the 
initial elastic response. We thus obtain surface strains increasingly 
concentrated near the rupture as time progresses. 

Observations of strains localized near the San Andreas 
fault following 1906 led Thatcher (1975a) to postulate his aseismic 
sliding model, but as shown here it is also true that such a 
strain "distribution can equally well be explained by anelastic 
adjustment. The important point is that whereas stable sliding , 
tends to alleviate overall strain buildup and thus increase 
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the time interval between shocks, viscoelastic adjustment has 
the opposite effect of increasing the stress on the fault following 
the earthquake and thus decreasing the waiting time until the 
next shock. 

As another- illustration, angle changes for the model tri- 
angulation networks in figure 10 have been computed from equa- 
tions (8) , (9) and (14) and the results are shown in Table 1. 

The most important factor determining transient motion due to 
viscoelastic adjustment within each time interval At = t = n/h 0 
is the relationship of the elastic layer thickness H to the 
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fault plane height D-d. If H'iD-d, the amount of post-rupture 
movement will be large; if H>> D-d, the motion will be small. 



Accuracy of the Green's Function 


To determine the accuracy of g ^ > we can compare the sur- 
face displacements computed from the approximate Green's function, 
equation (6) , with those obtained from a numerical integration 
of the exact solution (Ben-Menahem & Singh 1968) . Note the 
error in equation (11-57) of Ben-Menahem and Singh (1968) r which 
is corrected in Jovannovich et al. (1974a) . The integration 
was performed using the IBM-supplied Scientific Subroutine Package. 

In presenting the results of the surface displacement cal- 
/ 

culations we adopt the Jovannovich et al. (1974b) convention 

2 

of U Q = 1 km and dS — 1 km ' . The worst agreement between equa- 
tion (6) and the exact solution is expected when the source 
point is near the layer/half-space boundary, because the dis- 
placements and stresses implicit in the approximate Green's 
function are not continuous there. Tables 2-3 and figures 
11 - 12 are examples of comparisons between exact and approxi- 
mate Green's functions. 

In figure 11 and table 2, the layer thickness has been 
chosen equal to 10 km. Using fault parameters appropriate to 
the San Andreas as in Bundle and Jackson (1976b) , use was made 
of Thatcher-' s (19 74b) result that the fracture of the 1906 San 
Francisco earthquake extended to a depth of 10 km. The elas- 
tic layer was therefore constrained to be at least 10 km thick. 
Since the fault plane extends from the surface to 10 km depth, 
it is natural to check the approximate Green's function for 
H = 10 km, since as stated above, the approximation is expected 
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to be worse for tins case. 

In figure 11 and table 2, the source depth is varied from 
10 to 1 km xn 1 km increments. The point of observation was 
located at a radial distance of 10 km from the epicenter of 
the source. The error is defined as { (approximate surface dis- 
placement-exact surface displacement) / (approximate surface 
displacement) }xi00%. It is evident that the greatest errors 
in the 1 approximate Green’s function occur when the source is 
near the layer/half-space boundary. Above and below a source 
depth of about 5 km the dxfference between exact and approxi- 
mate Green’s functions has opposite sign. Upon integrating 
the Green's function over the entire 10 km depth, the errors 
tend to compensate in this case. 

To estimate the error in the displacements from a finite, 
-rectangular fault obtained using equation (6) , we can sum the 

entries in column 2 of table 2 and divide them by the sum of 

the entries in column 3. The result of this calculation indi- 
cates that the discrepancy xn the displacements caused by using 
(8) , (9) 'and (14) xn place of the exact solution xs small. 

Figure 12 and table 3 illustrate the effect of a changing 
elastic layer thickness H upon the surface displacements due 
to a source fixed at 10 km depth. The point of observation 
was again loeared at an epicentral distance of 10 km, while 

H was varied from 10 to 30 km in 5 km increments. As shown, - 

the error in the surface displacements decreases sharply as 
the interface is removed to a greater depth. For thicknesses 
H greater than 15 km, the error is never greater than 5% in 



magnitude . 


Note that equation (6) xs essentially a perturbation expan- 

li — 

sion in the quantxty - — e £ about the homogeneous half-space 
Green's function for y 3 <H, x^<H, and xs a sxmilar expansion 
multiplxed by the constant 2y 2 /{y 1 + y 2 ) for x 3 <H, Y 3 >_H. From 
the observation that successive terms represent more deeply 
burred "sources" of the same strength/ thus giving rise to smaller 
contributions to the surface displacements/ we can write the 
series in the form 


h = F 1 ^' 5jHi+ z r X ) + b <*, X ) 

m=l i— 

1 00 „ [— " 

9 12 = „ + X ){1+ h 5 lV«' # 1 

j- « in — * _L — — 


} Y 3 <H, x 3 <H 

(16) 


y 3 >H, x 3 <H 


where F^(;g, y) is the Green's function for a source in a homo- 
geneous half -space 7 and A and B m are defined as in (6) . Note 
that y) and B^(y, y) are positive definite since they 

are derived from terms under the summation in (6) which all 
have the same sign as F^^c, y) • A m (y, y ^ and ^ are also 

bounded, and because j?| < 1 for all physical earth models/ both 
series in (16) converge (e.g. , Whittaker & Watson 1927, o. 11). 

We can show that the representation (6) has the correct 
Ixmiting properties. For the case = u 2 r q\ 2 and g^ 2 reduce 
to the Green's function for the homogeneous half-space. Also, 
we know that integratxon of the solution for a point source 
along a line gives rise to the solution for a line source, so 
integration of the infinite series of point sources represen- 
ted by (6) along a line parallel to the x x axis yields Rybicki's 
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(1971) Green’s function. Finally , as we let H go to infinity 
in the upper of (6) , the terms under the summation vanish, whereas 
if H approaches zero in the lower of (6) we get 


lim g:L(x, 
H^O x 




2|i, 


y l~^2, 


m 


x ){1 + t ) 


+ V- 


2u- 

u7tn; F i ( ^ ^ 


m=l 

1 


(17) 




= V?' V 


1 - 


tl+U 2 


The fact that the solution to this problem can be repre- 
sented fairly accurately as a rapidly converging series in the 
parameter 5 suggests a method for the evaluation of exact inte- 
gral solutions for problems of this type. At present, the inte- 
gration of the exact Green’s function has been done numerically 
by approximating the denominator as a finite series of exponen- 
tials (Ben-Menahem & Gil Ion 1970) . Upon examination of the 
exact Green's function integral, equation 11-56 of Ben-Mehanem 
and Singh (1968) , we see' that if we write 
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£ = 


v 


1 - v 

1 + V 


1+5 


W HS = v ^] 


J HS 


X + y v 

X + 2y T 


x + u l ; rr| : 
x+ 3 Vtt|> 


( 18 ) 


we can expand the integrand in powers of 5* The convergence 
properties of the series can be investigated and if (6) is an 
indication, term-by-term integration can probably be carried 
out. The advantage of this procedure is that it may well put 
the integral in a form which is more amenable to analytic solution. 
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We know in advance that the first term must be the Green’s function 
for the homogeneous half-space and thus we know its algebraic 
representation (Steketee, 1958) . In addition , the modelling 
of an elastic layer over a viscoelastic half-space via. the cor- 
respondence principle is particularly easy due to the isolation 
of the elastic constants in the factor £ . One can extend this 
procedure to half-spaces of many layers by using multiple sums 

y i 

of several parameters £. = — where y is the rigidity of a 

y o ° 

particular reference layer. 



Summary 


V7e constructed an accurate analytic approximation to the 
Green’s function for a vertical strike-slip fault m an elastic 
layer over an elastic half-space. To model the San Andreas 
fault, we integrated the Green's function over a rectangular 
surface and introduced linear viscoelasticity into the half- 
space by means of the correspondence principle. Examples of 
viscoelastic relaxation in the form of post-fracture surface 
displacements and angle changes from a model triangulation net 
were computed. 
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Tr- 
iable 1. Comparison of the coseismic elastic and postseismic 
viscoelastic angle changes for two different ruptures, each 
in an elastic layer overlying a viscoelastic half-space {see • 
figure 10). For case I, the layer thickness is 20 km, the fault 
length is 400 km, the displacement dislocation across the fault 
face is 4 m and the fault has fractured the entire thickness 
of the layer. For case II, the layer thickness is 5 km, the 
fault length is 400 km, the displacement dislocation across 
the fault face is 4 m and the fracture extends from the surface 
to 3 km depth. 


CASE I 



CASE II 


Angle 

Coseismic 

Postseismic 
viscoelastic 
change 
(t = 2 t) 

Coseismic 

Postseismic 
viscoelastic 
change 
<t - 2x ) 

#' 

(sec) 

(sec) 

(sec) 

(sec) 

1 

■-6.16 

-1.83 

-3.47 

-1.10 

2 

-13.45 

-.98 ■ 

-9 . 84 

-1.12 

3 - 

9.59 

2.60 

6.40 

1.44 

4 

- 3.55 

-2.03 

5.58 

-1.97 

$ 

3-87 

-1.71 

4 ! .50 

-1.37 

6 

-7-43 

/ 

-3.75 

-10.09 

3.35 

7 

-.33 

-.04 

-.01 - 

-.02 

8 


-.13 

.67 

-.19 

9 

-.0 T 

.18 

-.66 

.22 

20l 

-12.32 

4.60 

- ^ 47 

10.55 

11 

la. 09- 

-5.85 

57.86 

-13.25 

12 

-3.77 

1.24 

-10.86 

2.70 
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Table 2. See figure 11. In column 2 are the horizontal sur- 
face displacements calculated from the approximate Green's function. 
Column 3 is the displacements computed with the approximate 
Greens ' function minus the displacements computed by a numeri- 
cal integration of the exact Green's function. 


Source 

Depth 

Approximate 

Displacements 

Displacement 

Difference 

Error 

(km) 

(cm) 

(cm) 

(%) 

10 

87.15 

12.36 

14.18 

9 

91.47 

11.99 

13.11 

8 

100.52 

9.13 

9.08 

7 

114.20 

5.60 

4.90 

6 

132.20 

.60 

.45 

5 

153.52 

-3.49 

-2.27 

4 

176.28 

-7.37 

-4.18 

3 

197.35 

-11.01 

-5.58 

2 

212.66 

-15.01 

-7.06 

~1 

”217.80 

-19.09 

-8.77 

Totals 

1483.15 

-16.29 

-1.1 
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Table 3. See figure 12. Columns 2-4 are the same as in table 2. 


Layer 

Thickness 

Approximate 
Di sp 1 ac eme n t 

Displacement 

Difference 

Error 

(km) 

(cm) 

(cm) 

(%) 

10 

87.15 

12.47 

14.31 

15 

52.04 

-.75 

-1.44 

20 

48.00 

-2.37 

-4.94 

25 

47.14 

-1.72 

-3.65 

30 

46.85 

-1.18 

-2.52 
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Figure 1. 
Figure 2 . 


Figure 3. 


Figure 4 . 
Figures 5 


Figures 8 


Figure Captions 


Coordinate system used in finding the displacements 
due to a line dipole of force in the y ^ direction. 
Coordinate system used to find the displacements 
due to a line dipole of force in the y ^ direction 
for a medium consisting of an elastic layer over- 
lying an elastic half-space. 

Coordinate system for the integration of the Green’s 
function for a strike slip point source. The rela- 
tive displacement across the rectangular fault sur- 
fact shown is indicated by the two parallel half- 
arrows. After Chinnery, M. , Bull, seism. Soc. Am. , 

^5, p. 9,, 1963/ copyrighted by the Seismological 
Society of America. 

A Maxwell element. 

7. Surface displacements parallel to a rectangular 
fault in an elastic layer over a viscoelastic half- 
space calculated using the approximate Green’s function. 
The displacements were calculated along the hori- 
zontal profile shown at the top of each figure where 
the double arrows lie along the fault trace. The 
Xayer thickness is 20 km and the fault fractured 
the entire thickness of the later. The relative 
displacement along the fault is 4m. 

9. Surface displacements perpendicular to a rectangular 
fault in an elastic layer over a viscoelastic half-space 



33 


All parameters are the same as for figures 5 - 7 . 

Figure 10- Location of some illustrative model survey triangles 
relative to a rectangular fault in an elastic layer 
over a viscoelastic half-space (see Table 1) • The 
fault is indicated by the double arrows, and the 
insert is a blow-up of the cross-hatched region. 

All angle changes were computed using the approxi- 
mate Green's function. 

Figure 11. Plot of the differences in the predicted horizon- 
tal surface displacements due to using the approxi- 
mate Green's function in place of a numerical integration 
of the exact Green's function for a point source 
in an elastic layer over an elastic half-space. 

H is fixed at 10 km and the rigidity of the layer 
is one-tenth the rigidity of the half-space. The 
point of observation is at a radial epicentral dis- 
tance of 10 km. See Table 2. 

Figure 12. Plot or the differences in the predicted horizon- 
tal surface displacements due to using the approxi- 
mate Green’s function in place of a numerical integration 
of the exact Green's function for a point source 
in an elastic layer over an elastic half-space. 

The rigidity of the layer is one-tenth the rigidity 
of the half-space and the point of observation is 
at a radial epicentral distance of 10 km. The source 
depth is fixed at 10 km. See Table 3. 





x Source Point at X3 
® Image Source Points 


Image Boundary 

II II 

y 2 

True Boundary 
Image Boundary 





PARALLEL OIS 



0 50 100 

PERPENDICULAR DISTANCE (km) 




PARALLEL DISPLACEMENT (cm) 


Fault 



0 50 100 
PERPENDICULAR DISTANCE (km) 



PERPENDICULAR 
DISPLACEMENT (CM) 


FAULT 

■ — - £ > 

1 5= p 

DISPLACEMENT.. PROFILE 

50 
40 
30 
20 

10 

0 10 20 30 40 50 60 70 80 90 100 

PERPENDICULAR DISTANCE (KM) 




PERPENDICULAR 
DISPLACEMENT (CM) 


FAULT 



DISPLACEMENT PROFILE 



0 10 20 30 40 50 60 70 80 90 100 

PERPENDICULAR DISTANCE (KM) 
















